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SUMMARY 

A general vector differential equation for the vorticity component 
parallel to a streamline is derived for steady, nonviscous and incom- 
pressible flow in a rotating system. This equation is then simplified 
by restricting it to rotating radial channels and by making further 
simplifying assumptions. This simplified equation is used to solve for 
the secondary vorticity, the vorticity component parallel to the stream- 
line, in three special cases involving different streamtube geometries; 
the results are presented in a series of figures. The secondary vorticity 
is shown to decrease with decreased absolute angular velocity of the 
fluid, decreased inlet total -pres sure gradient, decreased length of rela- 
tive flow path, and increased relative velocity. 


INTRODUCTION 

When a fluid with a total -pressure gradient normal to the osculating 
plane of the streamline follows a curved path, velocities are induced in 
the surfaces normal to the irrotational or "through" flow direction. 

These induced velocities or secondary flows are recognized as the source 
of several types of loss in turbomachinery. For example, losses may 
occur because of (l) boundary -layer separation caused by the transfer of 
low-energy air to regions of decelerating flow, (2) improper angles of 
attack in the blade rows of compressors and turbines, and ( 3 ) viscous 
dissipation and unrecoverability of the kinetic energy involved in the 
secondary flow. 

Because secondary flows are an important source of loss, considerable 
analytic work has been done on this phenomenon in stationary curved chan- 
nels by Squire and Winter (ref. l), Hawthorne (ref. 2), Kronauer (ref. 3), 
and others. The general method used in these analyses is to compute the 
vorticity component parallel to the actual streamline, the secondary vor- 
ticity, caused by a variation in inlet total pressure across the passage 
normal to the osculating plane of the streamline. This secondary vor- 
ticity is computed because it is an indication of the magnitude of the 
secondary flow. In addition, Kronauer (ref. 3) using a similar method 
has conputed this vorticity component in an axial-flow rotating channel 
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such as that in an axial -flow compressor. However, no work has been done 
on radial-flow rotating channels, such as in a centrifugal- type impeller. 

This report extends the work of Hawthorne and the others to a rotat- 
ing radial channel. The results of the analysis are indications of the 
qualitative trends of variables which affect the secondary vorticity and 
therefore the secondary flow. The purpose thus is not to obtain exact 
data but rather to obtain some insight into the problem of secondary flows 
in rotating radial channels. It is hoped that determining the relative 
importance of some of the variables involved will lead to a better under- 
standing of the flow in centrifugal-type impellers and similar flow ma- 
chinery. Calculations are carried out to show the effects of these vari- 
ables in several streamtubes with various geometries. 


METHOD OF ANALYSIS 

In this analysis the distribution along a streamline of the secondary 
vorticity (component of the vorticity parallel to the streamline) is com- 
puted for rotating radial channels. 


Preliminary Considerations 

For steady three-dimensional flow of a nonviscous fluid with nonuni- 
form total-pressure distribution through a stationary channel, the stream- 
lines and vortex lines lie on surfaces of constant total pressure (ref. 4, 
p. 244) called Bernoulli surfaces. The total -pressure gradient VP is 
necessarily normal to the Bernoulli surfaces, and in general these sur- 
faces are curved in space. The vorticity vector £, which is tangent to 
the vortex line, is defined by 

T = vxv 

where V is the absolute velocity vector tangent to the stre aml ine. From 
the equation of motion the vectors V, f, and VP are re la ted by 

vxr = vp 

Since VP is normal to the Bernoulli surface, this equation shows that 
the streamlines and vortex lines lie on Bernoulli surfaces. In general, 
the vorticity vector £ has components tangent and normal to the absolute 
velocity vector V, £ cos i|r and £ sin \|r, respectively, where \|r is the 
angle between the vortex line_and the streamline (see fig. l). The vor- 
ticity component tangent to V is the secondary vorticity, and it is this 
component of vorticity that is an indication of the magnitude of secondary 
flow. Between any two points on the Bernoulli surface there exists a 
geodesic or path of minimum distance. In general, the streamline is not 
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a path of minimum distance on the Bernoulli surface, that is, the stream- 
line is not a geodesic. If the streamline is not a geodesic, the princi- 
pal normal N of the streamline is inclined at an angle cp (fig* 1) 
from the normal to the Bernoulli surface VP* 

If the streamline is not a geodesic on the Bernoulli surface, the 
secondary vorticity increases in the direction of the flow (ref. 2). 

This condition exists if cp is not zero, that is, if VP has a component 
normal to the osculating plane of the streamline, the plane containing 
the velocity vector V and the principal normal N. If the streamline 
is a geodesic, VP lies in the osculating plane and the secondary vortic- 
ity does not change along the streamline. 

The general assumptions and method used in the development of the 
equations for secondary vorticity in a rotating radial channel are out- 
lined in the following sections. 

General assumptions . - The flow is assumed to he steady, incompres- 
sible, and nonviscous. Although the fluid is assumed to he nonviscous in 
the region of the solution, a gradient in inlet total pressure normal to 
the flow planes is assumed that would have been caused by the generation 
of a thick boundary layer by viscous forces upstream of the region of so- 
lution, that is, upstream of the rotating channel. If it is assumed, as 
will be done herein, that the boundary-layer thickness is large compared 
with the passage height (i.e., viscous shear forces are low), then the 
fluid can be assumed nonviscous and the vorticity calculated from the 
equation of motion for an ideal fluid. This assumption is also made in 
references 1 to 3. 

The flow entering the passage is assumed to be parallel and uniform 
except for the thick boundary layer which results in a total -pressure 
gradient normal to the osculating plane of the relative streamline. In a 
centrifugal-type inpeller, this would be the boundary layer on the inner 
and outer shrouds. Thus, initially, the vorticity conponent parallel to 
the streamline is zero. A conponent of vorticity parallel to the stream- 
line is generated when the streamline turns from its initial direction or 
is rotated in a rotating system. 

Outline of method . - In order to evaluate the secondary vorticity 
along a streamline, the vorticity component parallel to a streamline is 
expressed as a function of the relative total vorticity and the relative 
velocity. This expression is simplified using the fact that the diver- 
gence of the total vorticity is zero. The total relative vorticity is 
eliminated from the equation by means of the equation of motion. The 
equation of continuity is used to express the gradient along a streamline 
of the vorticity conponent parallel to a streamline divided by the rela- 
tive velocity as a function of the relative velocity, the angular velocity 
of the rotating system, and the gradient of the Bernoulli constant in the 



4 


MCA TN 3013 


rotating system. This is the general vector differential equation for 
secondary vorticity in a rotating system. This general vector differential 
equation is then simplified by restricting it to rotating radial channels 
and by making further simplifying assumptions. In order to solve for the 
secondary vorticity by means of this equation, it is necessary to know the 
streamline pattern. As the actual streamline configuration is, of course, 
unknown at the outset of the problem, it is necessary to assume a realis- 
tic streamline pattern which from general considerations would probably 
not differ greatly from the actual pattern in order to obtain meaningful 
results. This assumed streamline configuration is based on a potential 
flow solution or some other approximation, and a first approximation to 
the secondary vorticity is obtained by means of the simplified differen- 
tial equation for secondary vorticity. The secondary vorticity is solved 
for in several rotating radial streamtubes with various geometries using 
the simplified differential equation. 


General Vector Differential Equation for Secondary Vorticity 

The various steps outlined in the development of the general vector 
differential equation for secondary vorticity are now presented in de- 
tail. 


Vorticity vector . - The total relative vorticity is first expressed 
as the sum of its components parallel and normal to the streamline. Let 
W be the velocity vector relative to the rotating system and the 
relative vorticity vector equal to VXW. (All symbols are defined in 
appendix A and are expressed in terms of dimensionless ratios.) The 
dimensionless' relative velocity W is a ratio of the dimensional 
relative velocity to some characteristic speed of the c hann el cor^ . The 

angular velocity of the system is cn and r t is a reference distance 

from the axis of rotation to some point in the channel. The vorticity 
Q and the magnitude of its component | parallel to a streamline are 
dimensionless, having been divided by cs. 

_ _ It can be seen from the expansion of the triple vector product 
(WXH)XW that 


□ g-w - | (wxOjxw 

w-w w-w 

From the definitions of the dot and cross products it follows that the 
terms on the right are equal to the components of vorticity parallel and 
normal, respectively, to a streamline. Since £ is the magnitude of the 
component of the vorticity parallel to a streamline, the preceding equa- 
tion can be written 
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5= Iff + ( w *_ n > XH (i) 

w w*w 

where the absence of the bar indicates a scalar quantity. Because the 
divergence of the curl of any vector is zero, 0 can be eliminated from 
the left side of equation (l) as follows: 

V *0 = V-(Vxw) = 0 
so that equation (l) becomes 

vi w - v.ixinm (2) 

w ww 

In order to simplify equation (2) the equations of motion and continuity 
are introduced. 

Equation of motion . - The equation of motion can be written (ref. 5, 
eq. (157) 


W xfi = 


VH + 


2cdXW 

CD 


(3) 


where H is the Bernoulli constant relative to a rotating system, defined 
symbolically as 


, p/p + «£ . L 

(oa\(.)^ ^ ^ 


where p/p is the ratio of the static pressure to the fluid mass den- 
sity, made dimensionless by dividing by (cur^.) 2 , and R is the dimension- 
less distance in the radial direction expressed as a ratio of r^, as are 

all linear dimensions. The quantity H is analogous to the total pres- 
sure P in the stationary channel in that the relative streamlines lie 
on surfaces of constant H in the rotating case (fig. 2) just as the 
streamlines lay on surfaces of constant P in the stationary case 
(fig. l) . From equation (3) it can be seen that 


W X(WXft) = W X VH + 


ZHxtOTW) = - 2 £ _ 

“ CD 


2 — W (5) 

CD V ' 


By means of equation (5), the right side of equation (2) can be expressed 
as a function of W, co/cn, and VH. 
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Equation of continuity . - In order to solve for the component of the 
gradient of £/W parallel to the streamline , the equation of continuity 
is introduced. This equation states 

V.w = 0 

Hence , 

V-i W.= W-V i (6) 

and, from equation (5), 

y 1 XjWX a „ V .WXVH + g. v ( 7) 

W*W W 2 u*T 2 

General differential equation for secondary vorticity . - Equations 
( 6 ) and (7) are now substituted into equation (2), and the following 
differential equation results: 


w-v = v- w X — 

W v W 2 


+ 


w-v 


-2o>W 

(jdW 2 


( 8 ) 


The first term on the right of equation (8) represents the secondary 
vorticity growth caused by two effects: (l) the turning of the stream- 

line caused by the rotation of the system, and (2) the turning of the 
streamline caused by its curvature in the relative system. In order to 
clarify the physical significance of equation (8), the first term on the 
right is expanded as shown in appendix B to yield 


v .wx^ 2 vh“ + 

y2 y2 CO W 


(9) 


where Kg is the geodesic curvature of the streamline on the Bernoulli 
surface. Combining equations (8) and (9) results in the following general 
vector differential equation for secondary vorticity: 


W-V i = 4 : VH-- + 


2K, 


cNI 


W 


W £ 


03 


w 


+ W-V 


-2co-W 

ooW 2 


( 10 ) 


This vector differential equation relates the secondary vorticity £ at 
a point along the streamtube to the angular velocity of the system 
the relative velocity of the fluid W, the geodesic curvature of the 
streamtube K , and the gradient of the Bernoulli constant VH. 

o 
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Physical interpretation of general vector differential equation . - 
The first term on the right of equation (10) is influenced by the com- 
ponent of VH parallel to the axis of rotation, and the second term is 
influenced by the component of VH normal to the principal normal of 
the streamline. The third term on the right of equation (10) represents 
the variation of the vorticity component parallel to the streamline 
caused by the rotation of the system, irrespective of any variation of 
H. This type of secondary motion is discussed in reference 6 for an 
axial-flow-type channel. 

From the first two terms on the right of equation (10) it can be 
seen that in a rotating radial channel such as in a centrifugal impeller, 
in which the osculating planes of the relative streamlines are approxi- 
mately normal to the axis of rotation, the, only component of the gradient 
of H which affects the generation of secondary vorticity is that paral- 
lel to the axis of rotation. This corresponds to the component of the 
gradient H from the inner to the outer shroud. Thus, the boundary 
layer on the blade surfaces with the concomitant gradient of H normal 
to the blade surface does not affect the generation of secondary vor- 
ticity. 

The vorticity component parallel to the streamline induced by the 
rotation of the system is equal to -2cd cos a where a is the angle 
between the streamline and the direction of the axis of rotation (see 
fig. 2). The applications in this report are to rotating radial channels 
so that this term is zero because a is 90°. If the channel is sta- 
tionary, that is, a) is equal to zero, equation (10) becomes 



W W 


where it is understood that the variables are in dimensional form be- 
cause cd is equal to zero. This is the same equation as that obtained 
by Hawthorne in reference 2 for a stationary curved channel. 


Application of General Vector Differential Equation to Rotating 

Radial Channels 

In this section the general vector differential equation for second- 
ary vorticity developed in the preceding section is applied to rotating 
radial channels similar to those in centrifugal-type impellers, an exam- 
ple of which is shown in figure 3, and further simplifying assumptions 
are made. Numerical examples are solved involving streamtubes. with var- 
ious geometries which follow logarithmic spiral paths, including a 
straight radial path as a special case. A logarithmic spiral path is 
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similar to the streamlines in a conventional centrifugal impeller and 
an analytic expression is known for such a path. 

Further assumptions . - A right-handed orthogonal curvilinear coor- 
dinate system with coordinates u-^, Ug , and u-j is chosen. The unit 

vectors in the u^-, Ug-, and Uj-directions are i, j, and k, respec- 
tively. The reference point denoted by the subscript t in the general 
analysis shall refer to the tip of the rotating radial channel so that 
all linear dimensions are expressed as ratios of the tip radius and all 
velocities are expressed as ratios of the tip speed our^. 


The flow surfaces on which the relative streamlines lie are 
Bernoulli surfaces. These surfaces are parallel to the u-^ug-surfaces 

and are assumed to be planes normal to the axis of rotation (fig. 4). 
As a result of this assumption. 


a) = ook 


and 


_ 1 _ dH - 

u 3 ^ u 3 


( 11 ) 


( 12 ) 


where the U n for n = 1, 2, 3 are the square roots of the components 
of the metric tensor associated with the orthogonal curvilinear coordi- 
nate system employed and are given by 


U, 


n 



(13) 


in which V is expressed in the coordinates in which the u n are given. 
From this assumption it follows that the Bernoulli surfaces are planes 
coincident with the osculating plane of the relative streamline through- 


out the channel and cp is +90° or -90 c 


, , * dH 

" hen U7 35T 


is positive or 


negative, respectively. This assumption, which is usual in theo- 
retical studies of secondary flow, implies that the secondary flows 
are small perturbations on the primary flow. Hence, if the secondary 
flows are large, the results of this analysis are qualitative rather 
than quantitative. 


The orthogonal curvilinear coordinates in the u-^Ug -plane are se- 
lected so that lines of constant Ug are streamlines. Hence, 


W = Wi 


(14) 
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It is not necessary for the application of the method that the coordi- 
nate system he so chosen hut is merely a matter of convenience. As 
stated previously, it is necessary to assume a streamline pattern in 
order to solve for the secondary vorticity. The accuracy of the solu- 
tion for secondary vorticity is therefore dependent upon the accuracy 
of the streamline approximation. This approximation might he a two- 
dimensional potential flow solution, a mean flow path through the chan- 
nel, potential flow through the stationary channel, and so forth. 

Although the Bernoulli surfaces are assumed to he essentially paral- 
lel planes, the gradient of the Bernoulli constant is allowed to vary 
inversely as the assumed streamtuhe thickness ratio T in order to give 
a closer approximation to conditions in an actual passage. The stream- 
tuhe thickness ratio T is the ratio of the streamtuhe thickness in the 
uj-direction at any point to its value at the inlet of the rotating 

radial channel. Thus, 


_1_ BH _ 1 / l_ dH \ 
U 3 du 3 T \U 3 du 3 


(15) 


Simplified differential equation for secondary vorticity in a rotat- 
ing radial channel . - Under the assumptions and definitions given hy 
equations (11), (12), (14), and (15), the differential equation (10) for 
secondary vorticity in a rotating radial channel becomes 


_i_ a ft_\ = 

Ui Su^ \ W / 


w 3 t 


1_ 

u. 


dH 

^3 


(1 + W k) 


or along a streamline , 


( i) = _S_ (± |2-\ (1 

W W 2 t Ws 


+ w K) 


UldUi 

w 


(16) 


The term £/W is proportional to the secondary circulation, that is, the 
circulation associated with the secondary vorticity because £ is the 
secondary vorticity and l/w is proportional to the streamtuhe cross- 
sectional area normal to the flow and therefore normal to the secondary 
vorticity vector. It can he seen from equation (16) that the differen- 
tial change in secondary circulation parameter d(|/w) is directly pro- 


portional to 


1 ( 1 _ dH \ 
T \U 3 Su^ 


which is the Bernoulli constant gradient. 
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Also the differential change in the secondary circulation parameter is 
proportional to the sum (l + W k) and to U^u-l/w. In the sum 

(l + W k), the term Wn is the angular velocity about the axis of 
rotation of the fluid relative to the rotating system, and the term 1 
is the angular velocity of the system itself. Thus, the sum is equal to 
the absolute angular velocity of the fluid particle. Because the factor 
Updui/W is the differential element of time required for the fluid 
particle to pass through the differential channel length du-^, the 

U-.du-. 

product (1 + W K) — y — is the differential absolute angle through 

which the fluid particle is turned. In addition, d(|/w) varies inverse- 
ly with W 2 so that large values of d(£/w) will occur for small 
values of W • 

If the velocity W is constant, the secondary vorticity % varies 
linearly with the secondary circulation parameter £/W. If the flow is 
accelerating or decelerating for a given differential change in £/w, 
the differential change in ||| is greater or less, respectively, than 
that for constant W, as can be seen from the second term on the right 
in the following expansion of d(£/w): 



However, because the differential change in secondary circulation param- 
eter d(£/W) varies inversely with W (see eq. (16)), the tendency of 
the first term on the right of equation ( 17 ) is to produce an effect 
opposite to that of the second term. That is, accelerating and decel- 
erating flows tend to decrease and increase, respectively, |£j. Thus, 
there is a double effect involved so that it is not possible, in gen- 
eral, to predict the exact effect of accelerating or decelerating flow. 

Logarithmic -spiral coordinates . - A logarithmic spiral curvilinear 
coordinate system is specified in the relative flow planes and a recti- 
linear coordinate normal to the flow planes. Hence, the equations of 
the coordinates are 


u-l = In R + 6 tan p 
U2 = 9 - tan P In R 


( 18 ) 


u 3 = Z 


J 
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where R, 0, and Z are the usual cylindrical coordinates forming a 
right-handed system and 3 is the streamtube angle, that is, the angle 
between the radial and uq-directions, positive when in the direction 
indicated by positive oo according to the right-hand rule (see fig. 4). 
For logarithmic spiral streamtubes, 3 is constant along the streamtube. 

The curvature K of a logarithmic spiral streamtube is sin 3/R 
and the differential streamtube length is sec 3 dR. Hence, along a 
logarithmic spiral streamline equation ( 16) becomes 

d ^)- z -Wr (i\ i 1 + “ (19) 

The following equation results from integrating equation (19) along a 
relative streamline and noting that is equal to zero because the 

flow was assumed to be initially parallel and uniform except for a total- 
pressure gradient normal to the osculating plane of the streamline: 



where x is an arbitrary point along the streamline. This equation can 
be used to compute the secondary vorticity at any point x along the 
streamline from the known geometric and flow conditions. 

Case I. Constant t and WR . - In Case I, 


w t 

w = (21) 

and the streamtube geometry is similar to the channel in a centrifugal 
impeller with constant blade height and blade spacing. In this case the 
central angle subtended by the streamtube width Ad is constant. Be- 
cause T is constant, Bh/SZ is constant along the streamline. Inte- 
grating equation (20) using W as determined by equation (2l) yields 

[zfe < R ‘ - R i 4 ) + sl ” 13 - R l 2 >] < 22 > 

Equation (22) gives the secondary vorticity £ at any radius R for 
Case I. It can be seen that £ is directly proportional to (3H/dz)q 
and that the change in £ between any two radii is dependent on the 
value of Rj. 
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Several special cases can be considered with this particular geome- 
try. If Ri is equal to zero, equation (22) becomes 


e 




+ R sin P 


(23) 


If the strearntube walls are straight radial lines, P is zero and equa- 


tion (23) becomes 



Also, if the channel is stationary, a) 




2w t 2 


is zero and equation (23) becomes 

ban p 
W 


All the quantities in the above equation must be expressed in their 
dimensional form because 0 ) is equal to zero. 

The condition for £ to be zero at some radius R is given by 

r 2 + p.2 

p = 0 (24) 


Case II, 


(25) 

the channel in a centrifugal 
is constant) and constant flow 
area (blade height varies inversely with radius). From equations (15) 
and (25), Sh/SZ is directly proportional to the radius R. Integrating 
equation (20) using equation (25) for T yields 

5 /( I), = Vif 1 [s? < r2 - r i 2) + sin » < R - R i>] < 26 > 

Equation (26) is the expression for secondary vorticity at any point 
in the channel for Case II. The change in | between any two radii is 
inversely proportional to the value of R^. 


• -i 


2W+ 


+ sin 


as indicated by equation (22). 

Case II. Constant tR and W . - In 

R. 


T = 


1 

R 


and the strearntube geometry is similar to 
impeller with constant blade spacing (A 9 
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From equation (26) the secondary vorticity is zero at any point in 
the channel when 


1 + 


W sin 3 

R25 


= 0 


(27) 


It is interesting to note from equation (27) that the secondary vorticity 
is zero at a point when the absolute angular velocity (left side of eq. 
(27)) at the mean radius between the inlet and that point is zero. Pre- 
viously, it was noted that if W is constant , the differential change 
in secondary circulation at a point is zero if the absolute angular ve- 
locity of the fluid was zero at that point. 


Case III. Constant W ^ and T . - In Case III both the relative ve- 
locity W and the streamtube height ratio T are constants so that the 
streamtube geometry is similar to the channel in a centrifugal impeller 
with constant blade height and constant flow area (blade spacing A 6 
varies inversely with radius). Integrating equation (20) having noted 
that W and t are constants yields 

1 1 (l\ - ^f 1 [i? < E - R i> + » ln j q] ( 28 > 


Equation (28) expresses the variation in \ along a streamline for 
Case III. It can be seen that the change in £ between any two radii 
is independent of the value of R-^. The condition for zero | at a 
point is 


1 + 1 P In iL = 0 (29) 

(R - Ri) Ri 

The average absolute angular velocity between the inlet and radius R is 
given by 



W sin (3 

1 + Tr -"%) 


In 


R 


so that equation (28) indicates that the secondary vorticity is zero at 
a point when the average absolute angular velocity (left side of eq. 
(29)) is zero. This case corresponds closely in geometry to that con- 
sidered by Hawthorne and others in stationary elbows. The condition 
for zero \ in Case III, namely, that the average absolute angular ve- 
locity of the fluid is zero, is analogous to the condition of zero turn- 
ing angle for zero secondary vorticity in a stationary elbow. 
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Summary of equations. - A s ummar y of the geometric and flow condi- 
tions as well as the final equations for secondary vorticity and the 
condition for zero | in each case is presented in the following table: 



Variable 

Case I 

Case II 

Case III 

Specified conditions- 

A0 

Constant 

Constant 

(A0) t 

R 

T 

Constant 

R i 

R 

Constant 

Sh 

3z 

Constant 

r /a h\ 

Ri \SZ) ± 

Constant 

w 

w t 

R 

Constant 

Constant 

Consequences 
in flow 

‘/(I), 

Eq. (22) 

Eq. (26) 

Eq. (28) 

Condition 
for |=0 

Eq. (24) 

Eq. (27) 

Eq. (29) 


RESULTS AND DISCUSSION 

Standard conditions . - The final equations express the secondary 
vorticity | as a function of (Sh/Sz)^, W.J., R, R^_, and 3. As has 

already been noted, the secondary vorticity is directly proportional to 
(Sh/Sz)^. The effects of the other variables are presented in a series 

of figures in which the variation of £/(dH/5z)i with two variables is 
shown with the other variables given their values at standard conditions. 
These standard conditions are : 

W t 0.4 

0 0 

R ± 0.4 

R 1.0 

Case I . - The variation of the secondary vorticity parameter 
£/ ( SH/ dZ ) ^ with radius ratio R for various values of the inlet radius 
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Ri and the streamtube angle 3 is shown in figure 5 for Case I. In 
figure 5(a) the effect of the inlet radius Ri on the growth of the 
secondary vorticity in passing through the channel is shown. The vari- 
ation of the secondary vorticity parameter with inlet radius is small 
for this case, so the choice of the value of the inlet radius for the 
standard conditions will not affect the presentation of the results to 
any considerable extent. 

From the curve of \ / ( Bh/Bz) ^ against R for various values of 

streamtube angle (3 in figure 5(h), it can be noticed that the change 
in secondary vorticity at the tip radius for a given change in stream- 
tube angle decreases as the streamtube angle decreases. This phenomenon 
is caused by the fact that although the absolute angular velocity com- 
ponent decreases as 3 decreases from 40° to -40°, the length of path 
through the channel at first decreases as 3 changes from 40° to 0 but 
then increases as 3 varies from 0 to -40°. As discussed following 
equation (16), the change in secondary vorticity varies directly with 
both the absolute angular velocity and the path length so that when 3 
decreases from 0 the effects partially cancel each other. For 3 equal 
to -40°, the secondary vorticity parameter becomes negative over a range 
of R following R^ because the absolute angular velocity is initially 

negative. 

Figure 6(a) presents the variation in the secondary vorticity param- 
eter at the tip |^./(Sh/Sz)^ with tip velocity ratio for several 

values of 3 • The effect of W^. is pronounced, as would be expected 
from equation (16). For equal to 0.2 and 3 equal to -40°, 

is greater than for W t equal to 0.2 and 3 equal to -20° or 0, because 

the elongation of the flow path more than balances the effect of reduced 
absolute angular velocity of the fluid. 

Cases II and III . - The variation of £t/(SH/Bz)^ with for 

several values of 3 is shown in figure 6(b) for Case II. As in Case 
I, the values of £ t /(BH/Bz)i for Wt equal to 0.2 and 3 equal to 
-40° are greater than the corresponding values for 3 equal to 0 and 
-20° for the reason mentioned in case I. Secondary vorticities are, in 
general, also much greater in this case than in Case I because for the 
same value of the gradient of the Bernoulli constant Bh/BZ at the in- 
let, the value of Bh/BZ at any other point in the channel is greater 
in Case II than in Case I. Again the marked increase in secondary vor- 
ticity with decreasing W t can be noted. Figure 6(c) presents the same 

information for Case III as is shown in figures 6(a) and 6(b) for Cases 
I and II, respectively. Again the secondary vorticity parameter in- 
creases rapidly for. small values of . 
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This same trend is evidenced in the experimental results reported 
by Spannhake in reference 7 (pp. 145-159). The apparatus used for the 
experiment was an S-shaped tube that was rotated about an axis as shown 
in figure 7. The fluid, water, enters and leaves in the axial direction 
along the axis of rotation. A loss coefficient X based on the total- 
pressure drop from inlet to exit was computed. A range of ratios of 
circumferential velocity of the pipe to through flow velocity of the 
water, corresponding to l/W, from 0 to 30 was used. The loss coeffi- 
cient for flow through this apparatus shows the same trend, that is, 
rapid increase as l/W increases, as that of the secondary vorticity 
in the theoretical case. Since the experimental work is only somewhat 
similar to the theoretical work, for in a channel such as the rotating 
pipe the streamline configuration is much more complex, and since vis- 
cous dissipation, not considered in the theoretical work, is an import- 
ant factor, the only conclusion which can be drawn is that the trends 
are similar. 

Comparison of Cases I, II, and III . - A comparison of the variation 
of l t /(^H/dZ) i with W t is shown by the solid lines in figure 8 for 

the three cases. It is immediately apparent that the secondary vorticity 
for Case I is less for all values of W-t than it is for Cases II and 
III. If the secondary vorticity parameter | t /(Sn/5z) i were plotted 
against inlet velocity Wj_, the curves for Cases II and III would be the 
same since the velocity is constant throughout the streamtube. The curve 
for Case I is shown as the dashed line in figure 8 for which the abscissa 
is W i . It is seen that the curve for Case I lies above the curves for 

Cases II and III. As discussed following equation (16), the effect of 
decelerating flow is to increase the value of d(£/w), because of the 
consequent smaller values of W, and at the same time to decrease d£ 
for a given value of d(£/w). Thus it appears that of these two effects, 
that of increasing d( £/W) and consequently d| predominates because 
£ is greater for Case I than for Cases II and III when all three cases 
have the same inlet velocity, that is, when the' velocity is less for 
Case I than for Cases II and III at all points in the streamtube other 

than the inlet. When all three cases have the same velocity at the tip, 

as shown by the solid lines in figure 8, the velocity is greater for 

Case I than for Cases II and til at all points in the channel other than 

the tip. Thus d(|/w) is smaller and also the flow is decelerating so 
that, as can be seen from equation (17), these two effects compound, with 
the result that the secondary vorticity parameter is smaller for Case I 
than for Cases II and III. 

The secondary vorticities for Case II are larger than those for 
Cases I and III because if (cffl/dz)^ is the same for all three cases, 

Sh/SZ at any other point is greater for Case II than for Cases I and 
III. The secondary vorticity, being directly proportional to c>H/oZ, is 
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therefore greater for Case II. All three curves show the same pro- 
nounced increase in secondary vorticity as W-t decreases. 


SUMMARY OF RESULTS AND CONCLUSIONS 

A general vector differential equation for the vorticity component 
parallel to a streamline in a rotating system is derived. This equation 
indicates that at a given point along the streamtube the secondary vor- 
ticity ^ is a function of the angular velocity of the system 0), the 
relative velocity of the fluid W, the geodesic curvature of the stream- 
tube K g , and the gradient of the Bernoulli constant VH. This equation 

is then simplified by restricting it to rotating radial channels and by 
making further simplifying assumptions. This simplified equation was 
used to solve for the secondary vorticity in three special cases involv- 
ing different streamtube geometries, and the results are presented in a 
series of figures. These figures and the equations from which they were 
obtained indicate that : 

1. The differential change in secondary circulation parameter £/W 
at a point in a rotating radial channel is directly proportional to the 
absolute angular velocity of the fluid, the gradient cffl/dz of the 
Bernoulli constant, and the differential element of time required for 
the fluid particle to pass through the differential channel length. The 
differential change in |/w also varies inversely with the square of the 
relative velocity. 

2. The secondary vorticities are lowest for Case I with constant 
streamtube thickness ratio T and streamtube width Ad because of the 
higher fluid velocities upstream of the channel tip. 

3. The secondary vorticities are highest for Case II with constant 
TR (where R is radius ratio) and W because the decreasing streamtube 
thickness downstream of the channel inlet results in increased Sh/SZ.' 

4. The secondary vorticity increased rapidly in all cases as the 
relative velocity at the tip decreased because of the inverse variation 
of d£ with the relative velocity. 

5. The decrease in absolute angular velocity of the fluid caused 
by a decrease in the streamtube angle' P caused the secondary vortici- 
ties to decrease except where the resulting increased path length offset 
this effect. 

6. The loss coefficient measured in an experiment with a rotating 
S-shaped pipe showed the same trend, that is, rapid increase as the 
relative velocity ratio W decreased, as evidenced by the secondary 
vorticity parameter in the theoretical results. 

Leiws Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 31, 1953 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

H Bernoulli constant relative to rotating system, dimension- 

less, eq. (4) 

i unit vector in u^-direction 

j unit vector in ug -direction 

k unit vector in u^-direction 

N principal normal of streamline 

P total pressure of fluid in stationary channel 

p _ static pressure, dimensional 
R radial length, dimensionless 

r radial length, dimensional 

U ,U ,U square roots of components of metric tensor associated with 
a ° orthogonal curvilinear coordinate system, eq. (13) 

u l> u 2> u 3 orthogonal curvilinear coordinates, eq. (18) 

V fluid velocity in stationary channel 

W velocity relative to rotating system, dimensionless 

Z axial distance, dimensionless 

a angle between tangent to streamline and direction of axis of 

rotation 

P angle between radial line and streamline, positive in direc- 

tion indicated by positive o> 

£ total vorticity in stationary channel 

9 coordinate in right-handed cylindrical coordinate system 

A 0 central angle subtended by streamtube width 


K 


curvature of streamline 
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K geodesic curvature of streamline 

o 

X loss coefficient 

| vorticity component parallel to streamline, dimensionless 

p fluid mass density, dimensional 

T ratio of streamtube height in axial direction at any point to 

that at inlet 

cp angle between principal normal of streamline and normal to 

Bernoulli surface, positive when indicating positive rota- 
tion about streamline according to right-hand rule 

i|r angle between vortex line and streamline in stationary channel 

Q total relative vorticity, dimensionless 

o> angular velocity of system, dimensional 

y vector operator 

Subscripts : 

n index equal to 1, 2, 3, denoting conponents in u^-, u 2 -, 

u^-directions, respectively 

i inlet 

t characteristic point in channel which is tip in case of rotat- 

ing radial channel 

x arbitrary point along streamtube 

Superscript : 


vector quantity 
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APPENDIX B 


EXPANSION OF V 


WXVH 
W 2 


In order to simplify equation (8) and bring out its physical sig- 


WX VH 

nificance, the term V — is expanded as follows: 


W £ 


V 


.WXJjpI = .( WxVH ) + -L vH -ft 
W 2 W 4 W 2 


(Bl) 


This expression is further simplified by expressing the first term on 
the right as 

2 

- . (WXVH) = -4 [(w-vw)-(w XVH) + (wxHHwxvh)] 

w 4 vr 

Since the streamlines lie on surfaces of constant H (see fig. 2), VH 
is nor mal to W and the preceding equation becomes 

2 

- w . (w XVH) = zL r(W-VW) . (WX VH) + W 2 VH * ft] (B2) 

W 4 W 4 


Equations (l) and (3) and the fact that VH is normal to W are used 
to simplify the term VH • ft as follows: 


VH- ft = VH 


ft 


W + 




0) 


XW +VHX W 


W • W 


-2VH • “ 

CD 


(B3) 


Combining equations (Bl), (B2), and (B3) yields 


V- 


WXVH 


W 


(W* W)-(W XVH) + •% VH • ^ 
W 4 VP 03 


(B4) 


The first term on the right_of equation (B4) involves a triple scalar 
product in which the term (W • VW) represents the relative acceleration 
of the fluid particle. This acceleration has two components , one tan- 
gent to the streamline and the other along the principal normal equal 

to W 2 k, where K is the curvature of the streamline. Hence, the 
triple product is equal to -W 3 |VH|K sin cp, where cp is the angle mea- 
sured from the principal normal of the streamline to the normal to the 
Bernoulli surface (in the direction of the gradient of H), positive when 
indicating a positive rotation about the streamline according to the 
right-hand rule. The quantity K sin cp is the geodesic curvature K g 
of the streamline, so that this term is zero when the streamlines are 
geodesics on the Bernoulli surface. Hence, equation (Bl) becomes 



NACA TN 3013 


21 


„ WXVH 

W # i ■ — ■ 

W2 


VH • - + 

CD 



(9) 


Equation (9) is used in the simplification of equation (8) to clarify 
the physical significance of the various terms of equation (8). 


i 
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Figure 1. - Vector diagram for stationary channel. 
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Figure 4. - Flow planes for simplified flow in rotating radial channel. 
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Standard conditions unless otherwise noted. 


(b) Effect of 0. 


Figure 5. 
Case I. 


- Concluded. Variation of | /(||^ with R 
Standard conditions unless otherwise noted. 


for 
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for Cases I, IX, and III. Standard conditions unless 
otherwise noted. 
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